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Abstract. Naturalness of electroweak symmetry breaking in the next-to-minimal supersymmetric 
model points to scenarios in which the lightest CP-even Higgs boson has SM-like ZZh coupling and 
t O ' mass m/, ^100 GeV and decays partly via /i — > bb but dominantly into two CP-odd Higgs bosons, 

! h — > aa , where m a < 2m/, so that a — » t + t~ (or light quarks and gluons) decays are dominant. These 

preferred scenarios correlate well with the observed excess of events at LEP for rrif, ^100 GeV. 
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In the minimal supersymmetric standard model (MSSM) a generic supersymmetric 
(SUSY) spectrum that leads to natural electroweak symmetry breaking predicts a mass 
for the Higgs boson which is already ruled out by LEP data. To satisfy LEP limits the 
SUSY spectrum has to be either very heavy or very special^] However in the next-to- 
minimal supersymmetric model (NMSSM) the Higgs boson can be where it is predicted 
from natural EWSB with a generic SUSY spectrum yQ and there is even a hint in LEP 
data that this might be the case 0]. 
^3 \ The NMSSM is a very attractive model. It provides an elegant solution to the jl 

problem of the MSSM via the introduction of a singlet superfield S. For the simplest 
possible scale invariant form of the superpotential, 

XSH u H d + ^S i , (1) 

the scalar component of S naturally acquires a vacuum expectation value of the order 
of the SUSY breaking scale, giving rise to a value of }l of order the electroweak scale. 
Its particle content differs from the MSSM by the addition of one CP-even and one CP- 
odd state in the neutral Higgs sector (assuming CP conservation), and one additional 
neutralino. We follow the conventions of [5]. The associated trilinear soft terms are 

XA k SH u H d + ^A K S 3 . (2) 

The final two input parameters are tan/3 = h u /hd and /i e ff = Xs, where h u = (H u ), 
h d = (Hj) and s = (S). These, along with niz, can be viewed as determining the three 
SUSY breaking masses squared for H u , Hj and S (denoted m^, and m|) through 
the three minimization equations of the scalar potential. Thus, the Higgs sector of the 
NMSSM is described by the six parameters A , K , Ax , A K , tan/3 , /i e ff. 
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Scenarios that lead to such a special SUSY spectrum were recently found, see for example mixed 
anomaly-modulus mediation (]]] or gauge mediation with gauge messengers yj]- Both scenarios generate 
large mixing in the stop sector which maximizes the Higgs mass allowing all experimental limits to be 
satisfied with a fairly light SUSY spectrum. 
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FIGURE 1. Left: fine tuning vs. the Higgs mass for randomly generated NMSSM parameter choices 
with tan j3 = 10 and M3 (mz) = 300 GeV. Green circles are the points that satisfy limits on SUSY spectrum. 
Blue squares either have m/jj > 1 14.4 GeV or, if m/,, < 1 14.4 GeV, large branching fraction for h\ ^a\a\. 
Orange stars represent scenarios in which h\ — > a\a\ is large and m aj < 2m/,. Right: expected and observed 

95% CL limits on C 1 ^ = [gzzh/SzzhsM^ 1 ^ ~ * me shown vs. m/,. Also plotted are the predictions 
for NMSSM scenarios that appeared in our original (low density) scan 13] that give fine-tuning measure 
F < 25 and m ai < 2m/, and that are consistent with the preliminary LHWG analysis code. 



We can use naturalness of EWSB to predict the mass of the Higgs boson. Let us define 
the measure of fine tuning as 

d log my 

F = Max p F p = Maxp 



dlogp 

where the parameters p comprise the GUT-scale soft-SUSY-breaking parameters. We 
randomly generate NMSSM scenarios with fixed tan/3 = 10 and M^{mz) = 300 GeV 
and plot fine tuning vs. the Higgs mass, see Fig. Q] (left). Scenarios that satisfy the 
experimental constraints on superpartner masses are labeled by green circles. Scenarios 
labeled by blue squares have in addition > 114.4 GeV and/or large Br(/*i — > a\a\) 
(typically, the lightest CP-odd Higgs boson is mostly singlet). If m fl| > Imy the later 
scenarios are further constrained by Zh — > Zaa —> Zbbbb limits and those with full 
strength ZZh coupling are ruled out unless > 1 10 GeV. However, there are no limits 
on the Zh — > Zaa — > Zx + x~ T + T~ final state for m h > 87 GeV and so the scenarios 
with m ai < 2m/, satisfy all experimental constraints. These scenarios (orange stars in 
Fig. Q]) are among the scenarios with the very lowest fine tuning one can find given 
the current constraints on superpartners. Naturalness of EWSB in NMSSM clearly 
points to scenarios in which the Higgs boson with SM-like ZZh coupling and mass 
m/jj ~ 100 GeV decays partly via hi — > bb but dominantly into two CP-odd Higgs 
bosons, hi — » aiai, where m ai < 2m^ so that ai — > T + T~ (or light quarks and gluons) 
decays are dominant. Further, since the ZZh and WWh couplings are very SM-like, these 
scenarios give excellent agreement with precision electroweak data. 

Quite interestingly, the largest excess of events in the search for the SM Higgs boson 
at LEP is for a test Higgs mass of m/, ~ 100 GeV [6]. This excess cannot be interpreted 
as the SM Higgs boson (a SM Higgs boson with a mass of 100 GeV would lead to 
~ 10 times larger excess). However, in the orange scenarios discussed above, the Higgs 
boson has SM-like ZZh coupling, but highly reduced Br(h — > bb) ; such scenarios can 



naturally explain the observed excess. In Fig. Q] (right) we plot C^y limits and C^ f 
predictions for NMSSM scenarios that appeared in our original (low density) scan @] 
that give fine-tuning measure F < 25 and m ai < 2m,},. All these points are consistent 
with the experimental and theoretical constraints built into NMHDECAY as well as 
with limits from the preliminary LHWG full confidence level/likelihood analysis (we 
thank P. Bechtle for doing this analysis). The eight points with the lowest fine tuning, 
F < 10, cluster near ~ 98 4- 105 GeV (see also Fig. 3 of yfl). A detailed description 
of these eight points, including precise masses and branching ratios of the hi and a\, can 
be found in Ref. 0]. A significant fraction of the F < 10 points with low fine tuning are 
very consistent with the observed event excess. 

Finally, let us discuss conditions under which scenarios with m ai < 2m^ and large 
Br(/zi — > a\a\) occur in the NMSSM. In the NMSSM, one of the CP-odd Higgses 
is massless in the Peccei-Quinn symmetry limit, K — > 0, or in the R-symmetry limit, 
A K ,Ai — > OD. The masslessness of a\ in the limit A K ,Ax —> can be understood 
as a consequence of a global U(1)r symmetry of the superpotential under which the 
charge of 5 is half of the charge of H u H d . In the limit A K ,A^ — > it is also a symmetry 
of the scalar potential. This symmetry is spontaneously broken by the vevs of H u , Hj 
and S, resulting in a (massless) Nambu-Goldstone boson in the spectrum. Soft trilinear 
couplings explicitly break U(1)r and thus lift the mass of the a\. For small trilinear 
couplings, the mass of the lightest CP-odd Higgs boson is approximately given as: 



where cos 6a measures the doublet component of the lightest CP-odd Higgs mass eigen- 
state, a\ = cos OaAmssm + sin OaAs- In the limit of large tan /3 or \s\ ^> v, cos 9a can 
be approximated by cos A ^ 2v/stan/3, which indicates that in these limits the lightest 
mass eigenstate is mostly singlet. 

Naively, an arbitrarily small mass for the a\ is achievable provided small values of 
A K and Ax are generated by a SUSY breaking scenario. However, A K and Ax receive 
radiative corrections from gaugino masses (A K only at the two-loop level) and we should 
naturally expect Ax(m z ) ~ ^(M 2 (m z )) > 100 GeV and A K (m z ) > few GeV. Much 
smaller values would require cancellations between the values of Ax, A K coming from a 
particular SUSY breaking scenario and the contributions from the radiative corrections. 

Eq. © indicates that the contribution from Ax ^ 100 GeV to is highly 
suppressed if the lightest CP-odd Higgs is mostly singlet. We define a mea- 
sure of the tuning in Ax(mz) and A K (mz) necessary to achieve m Ql < 2m& as 
Fmax = max { | d log rrr a /d log Ax (mz) \ , (Jlogm^ jd log A K {mz) \ } which is to be 
evaluated for given choices of Ax(mz) and A K (niz) that yield a given m ai [8]. This defi- 
nition reflects the fact that m ai is completely determined by these parameters for fixed 
A , K", /i ef f , tan /3 . The dependence of Fmax onA K and Ax is given in Fig. [2l As expected, 
the smallest fine tuning or sensitivity is achieved for as small A K and Ax as possible. Of 
course, as we discussed earlier, very small values of A K and Ax would require cancella- 
tions between the bare values and the RGE-induced radiative corrections and this kind 
of cancellation would not be visible from the definition of Fmax - However, this is not a a 
concern given that very small values of A K and Ax do not in any case lead to sufficiently 
large Br(/zi — > a\a\) that Br(/zi — > bb) is adequately suppressed, see Fig.[2l Thus, in the 
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U smT =300 GeV, n=150 GeV, tan|8=10. m^>90 GeV, m Sl <2m^ 
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x = allowed by LEP 
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FIGURE 2. Tuning in vs. A K (left) and Br(/ii -> fliai) vs. (right) for jU = 150 GeV and 
tanj3 = 10. Light grey (cyan) large crosses are points that satisfy all experimental limits. The dark (blue) 
diamonds are those points that do not have large enough Bi(hi — > a\a{) so as to suppress Br(/;i — * bb) 
sufficiently to escape LEP limits. 



region of parameter space where soft trilinear couplings are at least of order the typical 
RGE-induced contributions, which is also the region where Br(/?i — ► a\a\) is large, 
the tuning of the soft trilinear couplings can be as small as &{5% — 10%). Finally, we 
would like to note that Fmax, typically substantially overestimate the magnitude of fine 
tuning with respect to GUT scale parameters. For scenarios in which gaugino masses 
and soft-trilinear couplings are correlated (not free parameters but rather calculable 
from a SUSY breaking scale) there is in principle no tuning associated with the SUSY 
breaking necessary to achieve m ai < Imj, 

In conclusion, we reemphasize that in the NMSSM the SM-like Higgs boson can be 
where it is predicted on the basis of natural EWSB and that the LEP event excess in the 
Z + b's channel for a reconstructed Higgs mass of ~ 100 GeV is consistent with a 
scenario in which the h decays mainly via /i — >• aa — > t + t~t + t~ (2m T < m a < 2m^) or 
4 jets (m a < 2m t ) leaving an appropriately reduced rate for h — > bb. We speculate that 
similar results will emerge in other supersymmetric models with a Higgs sector that is 
more complicated than that of the CP-conserving MSSM. This suggests that the strategy 
for Higgs searches at the Tevatron and LHC should be modified. 
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